Fine Grid Asteroseismology of G117-B15A and R548 

A. Bischoff-Kim, M.H. Montgomery, D.E. Winget 

The University of Texas at Austin, Astronomy Department, 1 University Station, CI4OO, 

Austin, TX 78712, USA 

agnesOastro . as . utexas . edu , mikemonOastro . as . utexas . edu , 
dewOastro . as . utexas . edu 

ABSTRACT 

We now have a good measurement of the coohng rate of G117-B15A. In the 
near future, we will have equally well determined cooling rates for other pulsating 
white dwarfs, including R548. The ability to measure their cooling rates offers us 
a unique way to study weakly interacting particles that would contribute to their 
cooling. Working toward that goal, we perform a careful asteroseismological 
analysis of G117-B15A and R548. We study them side by side because they 
have similar observed properties. We carry out a systematic, fine grid search for 
best fit models to the observed period spectra of those stars. We freely vary 4 
parameters: the effective temperature, the stellar mass, the helium layer mass, 
and the hydrogen layer mass. We identify and quantify a number of uncertainties 
associated with our models. Based on the results of that analysis and fits to the 
periods observed in R548 and G117-B15A, we clearly define the regions of the 4 
dimensional parameter space ocuppied by the best fit models. 

Subject headings: Dense Matter — Stars: Oscillations — Stars: Variables: Other 
— Stars: White Dwarfs 



1. Astrophysical context 



G117-B15A and R548 are pulsating white dwarfs with atmospheres dominated by hy- 
drogen. These stars are called DAVs or ZZ Ceti stars. Their effective temperatures are 
close to 12,000 K. They are non-radial, g-mode pulsators, where the restor ing force is buoy- 



ancy. Because of its potential use to study axions and other exotic physics (llsern et al.lll992 



Kepler et al.ll2005al ). G117-B15A was the object of several asteroseismological studies. R548 



is similar to G117-B15A both spect roscopically ari d pulsationally so studying them side by 
side may provide additional clues. iBradley Jl998h performed such a parallel study. More 
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recentlv iBenvenuto et al.l (|2002[) published an asteroseismological fit to G117-B15A's period 
spectrum. IBenvenuto et al.l evolved models from the Main Sequence and included detailed 
physics, such as time-dependent diffusion. Their best fit model has an effective temperature 
and a mass consistent with spectroscopi c results, but the periods of the model do not match 
the observed period spectrum as well as iBradleyl 's fits do. 



We take a new approach to white dwarf asteroseismology, enabled by the computing 
resources available today. Using models similar to iBradleyl 's. we perform a systematic, fine 
grid search for the best fits to the period spectra of G117-B15A and R548. We study both 
stars side by side. While we do not ignore the spectroscopy entirely, we take a naive approach 
and base much of our results on the period fitting alone. At the end we compare our results 
with the spectroscopy and previous studies. I n the sam e naive spirit, we treat the helium 
layer mass as another free parameter. iBradleyl (119981 ) and IBenvenuto et al.l (120021 ) both fixed 



its value at 10 ^M*, where is the total mass of the model. 

We organize our paper as follows. In section 2, we present the observables we have for 
G 117-B15A and R548. In secti on 3, we s umm arize concisely the asteroseismological results 
of IBenvenuto et al.l (120021 ) and iBradleyl (119981 ). In section 4, we explain our method. We 
include a discussion of our models and a quantitative study of what input parameters matter 
most. We prese nt our re s ults in sect ion 5, along with a disc ussion of how they compare with 
the fits done by IBradleyl (119981 ) and IBenvenuto et al.l (120021 ) . We summarize and discuss our 
results in section 6. 



The stars 



G117-B15A and R548 have much in common and as we shall see, studying them side 
by side helps constrain the model parameters. In Table [1] we list the spectroscopically 
determined effective t emperatures and gravities of both stars, along with the modes we 
observe in those stars (|Mukadamll2004l ). These are the clues we have at our disposal to find 
best fit models to G117-B15A and R548. 

The most solid observables listed in Table [1] are the periods. For both G117-B15A and 
R548, periods are determined to better than a second. The 215s mode in G117-B15A and 
the 213s in R548 are very stable and determined to extremely high accuracy. By comparing 
the amplitude of G117-B15A's 215s mode in G117-B15A in the ultra violet to those in the 

1 mode. 



Robinson et al.l (119951 ) were able to determine confidently that it was an 



visible, 

R548's higher amplitude modes are close to the 215s and 271s modes in G117-B15A and are 
likely I = \ modes as well. 
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Table 1: Observed properties of G117-B15A and R548 



G117-B15A R548 



TeflF 


11620 ± 200 (a) 


11990 ±200 (a) 




12830 ± 350 (b) 






12250 ± 125 (c) 






12010 ± 180 (d) 


11865 ± 170 (d) 




12375 ± 125 (e) 




logg 


7.97 ±0.05 (a) 






7.92 ±0.05 (b) 


7.97 ±0.05 (b) 




8.10 ±0.15 (c) 






7.94 ±0.17 (d) 


7.89 ±0.17 (d) 


Periods [s] ) 


215.22 (19.8) 


212.77 (4.1), 213.13 (6.7) 


(Amplitudes [mma]) 


270.86 (7.1) 


274.25 (4.1), 274.78 (2.9) 




304.15 (8.8) 








187.27 (1.0) 






318.08 (0.9) 






333.65 (1.0) 



References. — (a.) iBereeron et al.l (Il995bl) : (b) [ Bergeron et alJ (jl995al ): (c) iKoester et al.l (|l994l ): (d) 
Koester fc Holberj |200l[ ): (e) iRobinson et all jlQdK 



Effective temperature and gravity are more poorly determined than the periods. We 
give an overview of the current s p ectros copy in Table [TJ To obtain an effective temperature 
for G117-B 1 5A. iRobinson et al.l (119951 ) assumed a log g of 7.97, following work done by 
Daou et al.l Jl99of L 



3. Previous fits to G117-B15A and R548 



We summari ze the b e st fit models iBenvenuto et al.l (|2002| ) found for G117-B15A in 
Table [2] and those iBradleyl (119981 ) found in Table [31 Mh is the hydrogen layer mass and Mhs 
the helium layer mass. We define $ in equation [H In essence, it is the average difference 
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between the observed periods and the calculated periods. The lower $, the better the fit. 



Table 2: Best fits to G117-B15A's observed period spectrum according to lBenvenuto et al, 
(l2002h . 



Mode identification A; = 1, 2, 3 *fc = 2, 3, 4 



Teff [K] 

M,/M0 
Mh/M, 

MHe/M, 

$ [s] 



11400 
0.500 

10-6 
10-2 



11800 
0.525 

1.48 X 10-^ 

10-2 

6.5 



Note. — * Benvenuto et al. picked this fit because it matched the spectroscopy of Koester fc Allard ( 200Cll ) 
better. iBenvenuto et al.l do not give a list of calculated periods for the fc = 1,2, 3 fit, but they show and 
state that it is comparable in quality to the fc = 2, 3, 4 fit. 



Table 3: Best fits to G117-B15A's observed period spectrum according to iBradleyl (Il998l ) 



Mode identification 


= 1,2,3 


A; = 2,3,4 


Teff [K] 


12160 


11460 


10790 


12530 


11860 


11120 


M./Mq 


0.55 


0.60 


0.65 


0.55 


0.60 


0.65 


Mh/M, 


3 X 10-^ 


2 X 10-^ 


1 X 10-^ 


1.5 X 10-^ 


8 X 10-5 


1 X 10-5 


Mne/M, 




10-2 






10-2 




$ [s] 


0.5 


0.8 


1.3 


1.1 


1.3 


2.5 



Both authors find two families of solutions; one where the mode identification for G117- 
B15A's three obs erved periods is fc = 1, 2, and 3 and the other where it is fc = 2, 3, and 4 
{i = 1). iBradleyi notes that for his models, both classes of solutions fit the observed periods 
and spectroscopic temperature equally well. Both authors find that the fc = 1, 2, 3 farn ily of 
fits have thin h ydrogen layers and the fc = 2, 3, 4 family thick hydrogen layers. iBradleyl 's fits 



are hotter than IBenvenuto et ahl 's. Those are two trends to keep in mind. 



Finally, we summarize the best fit models 



Bradley 



found for R548 in Table H At the 
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time, iBradleyl had at his disposal only two confirmed modes (the high ani plitude doublets 
near 213s and 274.5s) and one likely mode (the 318s mode). Again. iBradleyl finds two classes 
of solutions, one with the three known modes identified as /c = 1, 2, and 3 and the other as 
k = 2,3, and 4 (all £ = 1). 



Table 4: Best fits to R548's observed period spectrum according to IBradleyl (119981 ) 



Mode identification 


A; = 1,2,3 


A; = 2,3,4 


Teff [K] 


12440 


11560 


10790 


12190 


11970 


11320 


MJMq 


0.55 


0.60 


0.65 


0.55 


0.60 


0.65 


Mh/M, 


3 X 10"^ 


1 X 10-^ 


5 X 10-* 


1.5 X 10-^ 


5 X 10-5 


5 X 10-6 


Mne/M, 




10-2 






10-2 




$ [s] 


3.2 


1.3 


0.3 


4.4 


4.6 


3.1 



In Tabl e H] we have re-evaluated $ based on the knowledge that the 318s mode is 
indeed real. IBradleyl could not rely on that mode and left it out. Including the 318s mode, 
we immediately learn something new: the k = 2,3,4 fits are significantly worse than the 
k = 1 , 2, 3 fits . If we do not include the 318s mode, the two classes of fits are equally 
good. IBradleyl concludes that the k = 2,3,4 fits are better because they agree with the 
spectroscopic temperatures better. He also computes the rotational splitting of the two 
modes and finds that they also appear to favor the k = 2,3,4 fits. 



4. Method 

In essence, the problem we have to solve in white dwarf asteroseismology is the simple 
minimization of a function (the average difference between the calculated periods and the 
observed periods) with n variables, where n is 4 in the present study. Those variables include 
TefT, M*, Mne and Mh. Expressed mathematically: 

<l>(Teff,M„MHe,MH) = (AP) = ^Sjli |p-i--p°b«| (1) 

where N is the number of observed periods and A is a normalizing factor to account for 
the fact that at higher effective temperatures and higher masses, the period spacing de- 
creases, artificially increasing our chances of finding a good period match. For the regions 
of parameter space considered in this work, A ~ 1. 
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The simplest way to minimize $ is to compute it for all conceivable values of the 4 
variables and pick the smallest value we find. But the number of times we have to evaluate 
$ can quickly become astronomical. To make matters worse, each evaluation of the function 
requires the full computation of a white dwarf model. The White Dwarf Evolution Code 
(WDEC), described in the next section, allows us to compute a large number of models in a 
small amount of time. But even with the WDEC, building a detailed map of $ over all the 
relevant parameter space has only recently become practical on a standard desktop machine. 



4.1. Models 



The WDEC evolves hot polytrope models from temperatures close to 100,000K down to 
the temperature of our choice. Models in the temperature range of interest for the present 
study are thermally relaxed solutions to the stellar structure equations. Each model we 
compute for the grid is the result of such an evolutionary sequence. 



The WDEC is described in detail in iLainb fc van HornI fll975[ ) and IWoodI (jl990|). We 
changed a few things since the work done by iBradleyl (119981 ) . We used smoother core com- 
position profiles and e xperimented with the more complex profiles that result from stellar 
evolution calculations (ISalaris et al.lll997l ) . We updated the envelo pe equation of s tate t ables 
from those calculated by 



Fontaine et al. 



(119771) to those given by 



Saumon et al.l (119951 ) and 



use OPAL opacities ( Iglesias fc Roger^ 19961 ). We also treated the abundance of element 



differently at the boundary between the helium layer and the hydrogen layer. We a s sume d 
equilibrium diffusion profiles, following the prescription given by lArcoragi fc Fontaine! (119801 ). 
but do not make the trace element approximation; this was sh own to yield results quite sim- 
ilar to those based on time- dependent diffusion calculations (lAlthaus et al.ll2003l ). We still 
treat diffusion at the carbon-helium transition zone as a free parameter. 

There are three parame ters asso c iated with the shapes of the oxygen (and carbon) core 
composition profiles used by IBradleyl (119981 ): the central oxygen abundance (Xq), the edge of 
the homogeneous carbon and oxygen core (qfm) , and the point where the oxygen abundance 
drops down to zero (qo). q is a mass coordinate, defined as g(r) = — log(l — M(r)/M^,), 
where M(r) is the mass enclosed in radius r and M^, is the stellar mass. We show an example 
of a basic oxygen abundance profile and those three parameters in figure [1] along with a 
Salaris-like profile. In figure [21 we show abundance profiles for a fiducial model (12,400K, 
0.62 Mq, Mhb = 10^^'^M*, and Mh = 10^ '''^M*), along with a model with a sharper carbon- 
helium transition zone. We shall use both models in the next section. From now on, the 
in the expression of the helium and hydrogen layer masses will be implicit. We adopted 
Salaris-like profiles for the carbon and oxygen abundances. 
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4.2. The Parameters that Matter Most 

In order to determine the relative importance of each factor, we computed periods for 
the fiducial model we described above (one of our best fit models for G117-B15A) and 
then varied one parameter at a time. Each time, we compared the resulting periods with the 
period spectrum of the fiducial model. In our analysis, we did not compare the full calculated 
period spectra, but instead picked the three periods that were closest to those found in the 
observed period spectrum of G117-B15A. This gives us an idea of what parameters matter 
most for fitting G117B15A and R548 in particular (though the same parameters in some 
order also matter for asteroseismological fits to other white dwarfs). 

We list the results in order of decreasing significance in Table |5l In column 1, we list 
the properties of the fiducial model, in column 2 the changes we made in each case, and in 
column 3 the average difference between the periods of the fiducial model and the modified 
model (equation [H with A = 1). We also detail the difference for each of the three periods of 
the fiducial model that match G117-B15A's observed period spectrum (216.3s, 270.9s, and 
304.6s). The last two lines in Table [5] refer to the MLT treatment of convection in the models. 
a is the ratio of the pr e ssure scale height to the m ixing length. ML2 ref ers to the prescription 
of Bohm fc Cassinelli ( 1971 ) and MLl to that of Bohm-Vitense ( 1958h. While convectio n is 



responsible for important non-linear effects in the light curves, (e.g. iMontgomeryl |2005| ) . it 
has very little effect on the periods of the modes excited in the models because it only affects 
the outer 10~^^ mass fraction of the models or less. 

Examination of how each change affects each mode individually reveals a few striking 
features. For instance, the 271s mode is strongly affected by a change in the core parameter 
Xo (AP271S = 7.18s), while the other two modes are not (AP216S = 0.58s and AP305S = 1.69s). 
This, and other features apparent in the last three co lumns of table 5 ca n be e xplained by 



examining the weight functions of the relevant modes. iMontgomery et al.l (120031 ) used them 
as a diagnostic of the effect of composition transition zones on mode trapping. In figure [3l 
we show the weight functions of the first three i = \ modes for the fiducial model of Table O 
In the second panel of figure [31 we see that the k = \ mode resonates strongly with the base 
of the He layer. This shows that it is extremely sensitive to the He layer mass, and very 
insensitive to the H layer mass, which is what we see from AP216S in Table [51 Looking at 
the third panel, we notice that the k = 2 mode predominantly samples the chemical profile 
in the core; its period change due to a change in C/0 profile is 5 times that of the other two 
modes. This explains why the 271 s mode is so sensitive to the core parameters. The last 
panels shows that the k = 3 mode samples both of these features, as well as farther out in 
the model near the base of the H layer. This mode is the most sensitive of the three to the 
H layer mass. 
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Table 5: Parameters that matter most for DAV models 





1 h Q Tl (Td 


/A P\ 


^-T 216s 


A Pr.^, 

^-T 271s 


305s 


= 0.62 M© 


= 0.56 Mq (-10%) 


16.7s 


+10. Is 


+ 18.1S 


+22.0S 


Salaris profiles 


Basic profiles 


15.2s 


-7.40s 


-32.3s 


-6.01s 


Tefr= 12400 K 


TefT= 11160 K (-10%) 


12.9s 


+12. Is 


+15.4S 


+11. 3s 


Mh = 10-^-^ 


Mh = 10-*^-^ (-10% in log) 


4.59s 


-0.273s 


+3. 29s 


+10.2S 


Mne = 10-2-3 


Mne = 10-21 (-10% in log) 


3.88s 


-5.19s 


+ 1.95S 


-4.50s 


Xo=0.76 


Xo=0.68 (-10%) 


3.15s 


+0.58S 


+7.18S 


+ 1.69S 


Soft C/He transition 


Sharper C/He transition 


2.61s 


-4.28s 


+0.713S 


-2.84s 


qfm=0.50 


qfm=0.45 (-10%) 


1.18s 


+0.0439S 


-2.70s 


-0.791s 


Saumon et al. (1995) 


Fontaine et al. (1977) 


1.08s 








envelope EOS 


envelope EOS 










a = 0.6 


a = 2 


0.0104s 








ML2 


MLl 


0.0000689s 









In order to keep our model grids manageable, we decided to vary 4 parameters. While 
Teff, M^,, Mh, and Mne have by far the largest effect on the periods, the parameters associated 
with the core abundance profiles and the shape of the carbon-helium transition zone (second 
block in Table [5]) can have a significant effect (as high as ~ 3 seconds) on the pulsation 
periods. Since we had to fix our core composition profiles, we decided to fix them to profiles 
predicted by stellar evolution (figure [2]). We also adopted the fiducial model's carbon-hehum 
transition profiles. 

In the asteroseismological fits detailed in section 5, we use the periods rounded to the 
second, and do not distinguish a model period that fit within 0.5 seconds from one that fits 
within 1 second. This allows us not to worry ourselves with the factors listed in the last 
block of Table [5] and other possible small effects. 
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4.3. Fine Grid Search 

We started by building a low resolution grid that covered a broad region of parameter 
space, guided by the spectroscopy (Table [1]). We varied masses between 0.46 M© and 0.80 
Mq and temperatures between 10800 K and 13000 K. We considered —3 < log((MHe) < —2 
and —8.4 < log(MH) < log(MHe) — 2. For that grid, we determined that step sizes of 200K in 
TefT, 0.02 M0 in M^,, and 0.4 in log(MHe) and log(MH) were sufficient to locate likely minima 
of $. Using the results of the broad grid search, we narrowed down our search to smaller 
areas of parameter space and built a finer grid. We found that a resolution of lOOK in Tcs, 
0.01 Mq in M,,, and 0.2 in log(MHe) and log(MH) was sufficient to clearly define the minima 
of $. The broad grid contains near 20000 models, the fine grid 6000 (it is more limited in 
parameter space). 

For G117-B15A, we further narrowed down the list of possible best fits by assuming 
that G117-B15A's modes were i = 1 modes, consistent with mode identification work done 



by lRobinson et al.l (Il995l ) and the fact that higher 1 modes suffer from geometric cancellation 
at the surface of the star and are likely to result in low-amplitude pulsations. For R548, we 
required that the two high amplitude modes, which also seem to be present in G117-B15A 
(see Tabled]), he i = 1 modes. We did not place any constraints on the identification of the 
other three modes. 



5. Results 

5.1. Mass and Effective Temperature 

We display the results of the fine grid search for G117-B15A and R548 side by side in 
the M=K — Teff plane in figure HI We explored all of the parameter space shown. The fine grid 
starts at 0.6 Meand above. We also indicate the spectroscopic temperature determinations 
for both stars, as well as masses derived from the gravities listed in Table [H using our models. 
Our best fit models appear systematically massive and/or hot compared to the spectroscopic 
results. We shall come back to this discrepancy in section 6. 

With our Is sensitivity cut-off due to small modeling uncertainties, we are unable to 
determine a unique point in parameter space that matches G117-B15A or R548. Instead, 
we find families of solutions. There is a tigh t correlation b etween the mass and the effective 



temperature, anticipated from earlier work (jBradleylll998l ). Decreasing either the tempera- 



ture or the mass decreases the spatial average of the Brunt- Vaisala frequency in a similar 
way, and therefore yields similar sets of periods. In general, our best fit models are more 
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massive and/or hotter than we would have expected from the spectroscopy. 

From their observed properties (Tabled]), we expected G117-B15A and R548 to be best 
fit by similar models. We do see that in our results. It also comes as no surprise that the 
£ = 1 identification requirement does not limit model fits to R548 as much as they do for 
G117-B15A. To obtain the dotted circles in figures HI we discarded fits that did not obey 
the £ = 1 mode identification for G117-B15A's 3 modes and for R548's 2 high amplitude 
modes (regardless of the quality of the fit). Because this first cut is based on the constraint 
of a fewer number of modes in the case of R548, it does not eliminate as many models. On 
the other hand, we have 3 additional modes to fit for R548 and they ultimately allow us to 
narrow down the best fit models to a small region of parameter space. We do not have that 
luxury with G117-B15A, as the 3 £ = 1 modes are all we have. We also obtain better fits to 
G117-B15A's 3 periods than we do to R548's 5 periods, as we would expect. 



5.2. Helium and Hydrogen Layer Mass 

For both stars, the helium layer mass appears to be fairly well constrained around 
4 X 10^'^. It is determined by the 215s mode in G117-B15A and the 213s mode in R548. If 
we do not include those modes in the fit, Mne is essentially unconstrained. For G117-B15A, 
if we change that mode by as little as 5 seconds (e.g. 215s to 210s), the helium layer mass 
changes from 4 x 10^^ to 5 x 10"^. If instead we leave out the 271s mode, the helium 
layer mass remains constrained. This result is readily understandable from figure [3] and the 
associated discussion of weight functions in section 4.2. 

In figure El we show the very best fit models ($ < Is for G117-B15A and $ < 1.5s 
for R548) in the — Tcfr plane and indicate their respective hydrogen layer masses. We 
discover that families of models with different hydrogen layer masses separate out in the 
M^. — Teff plane. 

We find that R548 is best fit with thin hydrogen layer models (Mh ~ 2 — 6 x 10^^). 
For G117-B15A, we find two well defined families, one at higher mass and lower effective 
temperature with Mh = 6.3 x 10"'' and one at lower mass and higher effective temperature 
with Mh 1 — 4 x 10~^. Previous iri vestigations have found both thick and thin solutions 



(jBradleylll998l : iBenvenuto et al.ll2002l ). Based on the incompleteness of our current under- 



standing of mass loss in the late stages of stellar evolution, we do no t beheve either set o f 



models can be ruled out by stellar evolution calculations. For instance, lAlthaus et al.l (120021 ) 
present the hydrogen abundance (Mh = 10""^) in DAVs as an upper hmit, as it would be 
reduced by the inclusion of mass loss episodes during the planetary phase, the extent of 



which is unknown. While none of our best fits truly have thick hydrogen layers, we shall 
refer to the two distinct families as "thick" and "thin" hydrogen layer fits. 



5.3. Mode Identification 



Previous asteroseismological studies of G117-B15A (jBradleylll998l : lBenvenuto et al. I l2002h 
identified the three modes as cons ecutive i = 1 modes, with the 215s mode being either a 
k = 1 or a = 2 mod e. iBradleyl found equally good fits with either mode identification. 
While iBenvenuto et al.r s best fit model was consistent with a /c = 1, 2, 3 mode identification, 
they also found good fit models with k = 2,3,4 modes. Among the 27 models in the fine 
grid that matched the observed periods to better than 1 second on the average, we found 
13 matched the observed periods with A; = 1, 2, 3 modes and 14 with A; = 2, 3, 4 modes. The 
former family of models all have Mh — 1 ~ 4 x 10 ~^ an d the latter Mh = 6.3 x 10"''. This 



is qualitatively in tune with what iBradleyl (119981 ) and IBenvenuto et al.l (120021 ) found (the 
/c = 1, 2, 3 fits have thin hydrogen layers and the k = 2,3,4 fits have thicker hydrogen layers). 

For R548, among the top 20 best fit models, we find that the most likely mode identi- 
fication (17 out of 20 models) is £ = 1, = 1, 2, 4 respectively for the 213s, 274.5, and 318s 
modes and i = 2, k = 4 and 8 for the 187s and 334s modes. Recall that we required the 213s 
and 274.5s modes to be £ = 1 modes, but did not place any constraints on the other modes. 
Two out of the three models that disagree with this mode identification are thick hydrogen 
models (while the 17 models with the most common mode identification are thin hydrogen 
models). The i = 2, k = 4 identification for the 187s mode is very robust (20/20). 

The main reason we are studying G117-B15A and R548 side by side is that they are 
observationally similar and we therefore expect them to also be structurally similar. If that 
is the case, then the mode identification results for R548 suggest that the correct mode 
identification for G117-B15A is k = 1,2 and 3, in favor of the thin hydrogen layer and lower 
stellar mass solutions (Mh ~ 1 — 4 x 10~®). 



6. Summary and conclusions 

We performed a systematic fine grid search for best fit models to the two DAVs G117- 
B15A and R548. We find best fit models for both stars between 11600 and 12700 K and 
between 0.59 and 0.66 Mq. In both cases, the region occupied by the best fit models in the 
M^, — TefT plane is related to the thickness of the hydrogen layer. Treating the helium layer 
mass as a free parameter, we discovered that the lowest period mode for each star (215s for 
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G117-B15A and 213s for R548) singlehandedly sets the helium layer mass to 4 X 10-3. Both 
stars are well fit with thin hydrogen layer models (10"^ < Mh < lO"®). For G117-B15A, 
we find a second family of solutions between 11300 and 12300 K and between 0.65 and 0.68 
Mq. Those fits have slightly thicker hydrogen layers (Mh = 6.3 x 10"''). 

For R548 we find a unique, robust mode identification. The dominant modes (213s and 
274.5s) are i = 1, k = 1,2 modes. The 318s mode is also an i = 1 mode, with k = 4. 
The last two modes (187s and 334s) are £ = 2, k = 4 and 8 respectively. For G117-B15A, 
we find two distinct families of best fit models. Models with Mh = 6.3 x 10"^ all have the 
same mode identification, namely A; = 2, 3, and 4 for the three observed periods (215s, 271s, 
and 304s). Models with Mh ~ 1 — 4 x 10"^ are consistent with G117-B15A's periods being 
k = 1,2, and 3. R548 and G117-B15A have similar observed properties and based on that 
fact, it is likely that they have similar structures. In this case, the second class of models 
appears better, but we cannot discard the first class of models based on that fact alone. 

By sampling parameter space systematically and homogeneously, we found that our 
models were offset in mass and temp erature compared to the spectroscopy (the models are 
hotter and/or more massive). iReidI (1996) measured the gravitational redshift of spectral 
lines for G117-B15A. The resulting mass is in tune with the spectroscopic mass (0.536 ± 
O.OIOM0). Both different equation of state tables and opacities could alter the structure of 
the models and possibl y modify our result s. We do not expect either to have a large effect. 



however. Just recently, ICassisi et al.l (120071 ) published new electron-conduction opacities that 



treat the par tially degenerate regime relevant in white dwarf envelopes better. According to 



Cassisi et al.l . the new opacities differ only by at most a factor of 2 in white dwarf envelopes 



and have very l ittle e ffect on the mass-radius relation. We also tried to revert back to the 



Fontaine et al.l (119771 ) envelope equations of state and found that this had a neghgible effect 



on the gravities in our models. 

Results from section 4.2 suggest a more promising avenue. They show that structure in 
the core and the shape of the helium composition profile at the carbon-helium interface can 
have a large effect on the pulsations periods of a model. We adopted core composition profiles 
from stellar evolution calculations and did not vary them. Our asteroseismological results 
suggest that perhaps we should try different core composition profiles. With precise enough 
determinations of G117-B15A and R548's mass and effective temperature from spectroscopy 
or other independent methods, one can turn the problem around and attempt to determine 
what composition profiles are needed in order to reconcile the asteroseismology with the 
spectroscopy. While preliminary results appear promising, a full investigation requires the 
clever analysis of a more extensive grid of models and we leave that for a future publication. 

Once we are satisfied that we have models that match both asteroseismologically and 
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spectroscopically, we can calculate P's for our models and use them together with the ob- 
served P's to constrain the energy loss rate due to any weakly interacting particles, such as 
axions. The fine grid approach allows us to formaly assess the uncertainties in our mod- 
els' parameters and to obtain tight constraints on the emission rate of weakly interacting 
particles in G117-B15A. 
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-log(l-Mr/M*) 



Fig. 1. — Basic oxygen abundance profile (solid curve) and the three associated parameters 
Xq, Xfm, and qo- The dashed-dotted curve is an example of a Salaris-like oxygen abundance 
profile. 
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Fig. 2. — Chemical composition profiles for the fiducial model (solid lines) and for a model 
with a sharper carbon-helium transition zone (dashed lines). 



-17- 



log(l-MyMj 
-5 -10 




$(0 



Fig. 3. — A plot of the weight functions of the first three £ = 1 overtones as a function of 
the coordinate $(r), the "normahzed buoyancy radius", for the fiducial model of Tabled 
We see that the k = 1 mode has its period mainly determined by the composition gradient 
at the base of the He layer, and that the k = 2 mode is most sensitive to the structure in 
the C/0 profile in the core. The k = 3 mode is sensitive to multiple features in the model 
and the only one sensitive to the location of the base of the hydrogen layer. 
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Fig. 4. — The valley of best fit models for G117-B15A and R548 in the M^, — Teff plane. 
The dashed circles mark the location of the subset of models that fit the ^ = 1 mode 
identification criterion (see text). Of those, the progressively filled-in circles indicate better 
and better fits ($ < 2.5s, 2s, 1.5s, Is respectively). The filled symbols with error bars 
indicate the spectroscopically determined temperatures and mass for G117-B15A and R548, 
according to the legend. We used our models to derive a mass from the gravities listed in 
Table [D 
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Fig. 5. — Best fit models in the M,,, — Teg plane for different hydrogen layer masses. For 
G117-B15A, we isolated the solutions for which $ < Is and for R548, those for which 
$ < 1.5s. The open circles correspond to "thick" hydrogen solutions (Mh = 6.3 x 10~^), 
while the filled circles correspond to thin hydrogen solutions (Mh ~ 1 — 4 x 10^^). 



